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Abstract 

We report on the abundance analysis of Li in solar-type (G-type main sequence) 
superflare stars which were found by the analysis of Kepler photometric data. Li is a 
key element to understand the evolution of the stellar convection zone which reflects 
the age of solar-type stars. We performed the high dispersion spectroscopy of solar- 
type superflare stars with Subaru/HDS, and conhrmed that 34 stars show no evidence 
of binarity in our previous study. In this study, we derived the Li abundances of these 
34 objects. We investigate correlations of Li abundance with stellar atmospheric 
parameters, rotational velocity, and superflare activities to understand the nature of 
superflare stars and the possibility of the nucleosynthesis of Li by superflares. We 
conhrm the large dispersion in the Li abundance, and the correlation with stellar 
parameters is not seen. As compared with the Li abundance in Hyades cluster which 
is younger than the Sun, it is suggested that half of the observed stars are younger 
than Hyades cluster. The measured value of nsini (projected rotational velocity) 
supports those objects are younger than the Sun. However, there are some objects 
which show the low Li abundance and slowly rotate on the basis of the estimated 
nsini and P (period of brightness variation). This result indicates that the superflare 
stars are not only young stars but also old stars like our Sun. In our observations, 
we could not hnd any evidence of Li productions by superflares. Further research on 
Li isotope abundances of superflare stars would clarify the Li production by stellar 
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flares. 
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1. Introduction 

Lithium (Li) is easily destroyed in the hotter region (p, a reactions ; ^Li : T > 2.5 x 
10® K, ®Li : T > 2.0 x 10® K) of stellar atmosphere. We can obtain the information about the 
evolution of convection layer from the behavior of Li abundance. The accurate spectroscopic 
determination of the Li abundance in young solar-type stars provide independent and reliable 
age diagnostics (e.g., Herbig 1965, Duncan 1981). In addition, accurate measurements of the 
Li abundance and, in particular, the ®Li/^Li isotopic ratio in stellar atmospheres are of crucial 
importance for addressing questions about the Big Bang nucleosynthesis (e.g.. Spite & Spite 
1982, Boesgaard & Steigman 1985), the chemical evolution of the Galaxy (e.g., Ryan et al. 
2001), mixing processes in stellar interiors, and the evolution of extra-solar planetary systems 
(e.g., Bouvier 2008, Israelian 2009, Gonzalez 2015). 

Solar flares are the most energetic explosions on the surface of the Sun, and are thought 
to occur by release of magnetic energy (Shibata & Magara 2011). Flares are also known to 
occur on various types of stars including solar-type stars. Among them, young stars, close 
binary stars, and dMe stars sometimes produce “superflares”, flares whose total energy is 10- 
10® times larger (10®^ —10^® erg) than the largest flares on the Sun (~ 10^^ erg) (Schaefer et al. 
2000). It has been recognized that the superflares will not occur on the slowly-rotating single 
solar-type stars. However, Maehara et al. (2012) analyzed the photometric data by the Kepler 
spacecraft, and discovered 365 superflare events on 148 solar-type stars that have the effective 
temperature of 5,100K < < 6,000K, and surface gravity of \ogg > 4.0. They also found 

that the brightness of superflare stars show the quasi-periodic variability with period of from 
one to a few tens of days. Those brightness variations can be explained by the rotation of a 
star with large starspots (Notsu et ah 2013b). This result clearly shows that we can obtain the 
accurate velocity of stellar rotation. It is important to measure the rotational velocity of stars 
which is thought to be a key parameter of stellar evolution. The rotation velocity of a star also 
reflects the age and activity of a star. Slowly rotating stars should be quiet and old stars like 
our Sun. 

As a subsequent study, Shibayama et al. (2013) found 1547 superflare events on 279 solar- 
type stars, and 44 superflares on 19 slowly rotating “Sun-like” stars (5,600K < Teg < 6,000K, 
log(? > 4.0, brightness variation period (P) > 10 days). Maehara et al. (2015) also found 187 

* Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory 
of Japan. 
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superflares on 23 solar-type stars using Kepler data with 1 min sampling. With these data, 
we studied the statistical properties of the occurrence rate of superflares, and found that the 
occurrence rate {dN/dE) of superflares versus flare energy (E) shows a power-law distribution 
with dN/dE oc E~°‘ where a ~ 2 (Maehara et al. 2012; Shibayama et al. 2013; Maehara et al. 
2015). It is interesting that this distribution is roughly similar to that for solar flares. These 
results indicate that it is quite likely that superflares would occur on the Sun. In addition to this 
research, Shibata et al. (2013) show the possibility of occurrence of superflares on the present 
Sun from the theoretical point of view. Notsu et al. (2013a) show the chromospheric activity 
of the solar-type superflare star (KIC 6934317) by high dispersion spectroscopy. They found 
high chromospheric activity from the core depth of Ca II and Ha lines of KIC 6934317. These 
results support the idea of existence of large starspots. They also found the low Li abundance 
in that star. Nogami et al. (2014) obtained spectroscopic data of the slowly-rotating superflare 
stars KIC 9766237 and KIC 9944137, and found that these superflare stars have very Sun-like 
atmospheric parameters including the Li abundance of the Sun. Wichmann et al. (2014) carried 
out high dispersion spectroscopy for superflare stars from the list of Maehara et al. (2012). They 
found that several of those stars are very young and active, but some other stars do not show 
the cause of occurrence of superflares. 

Notsu et al. (2015a,b; hereafter Paper I, II) observed a large sample of “solar-type” 
superflare stars. They picked up single stars and conhrmed the atmospheric parameters of 
“solar-type” superflare stars are similar to the solar ones by high dispersion spectroscopy. In 
particular, they found that 9 stars are in the range of “Sun-like” stars. This result supports 
the hypothesis that the Sun might cause a superflare. In paper II, they showed the correlation 
between the intensity of Ca II infrared triplet lines, which are good indicators of the stellar 
chromospheric activities, and amplitude of brightness variation among their observed superflare 
stars. These results clearly show that all the targets have large starspots and high chromo¬ 
spheric activity compared with the Sun. From those observational results, we can say that the 
brightness variation of snperflare stars is due to the rotation with large spots. The problem 
which we have to consider next is whether the present Snn causes the superflare or not. 

Several studies have been done on the behavior of Li abundances in solar-type stars so far. 
Takeda et al. (2007, 2010) observed the large sample of “solar-analog” stars (selected with the 
photometric criteria of0.62^i? —K ^0.67 and 4.5 ^ My ^5.1) and investigated the correlations 
of stellar parameters with Li abnndance. They fonnd clear correlations among Li abundance, 
stellar activity and rotation. Mishenina et al. (2012) also investigated the correlations among 
Li abnndance, stellar parameters and activity in the F,G,K dwarfs. They show the Li-activity 
correlation is evident only in a restricted temperature range (5,200 < < 5,700K), and the 

spread of Li abundance seems to be present in a group of low chromospheric activity stars. 

The possibilities of Li production in stellar flares have been suggested by some researches. 
Canal (1974) shows the possibility of nucleosynthesis of Li in low-energy flares in T Tauri-like 
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stars. Karpen & Worden (1979) estimated the amounts of Li produced by stellar-flare processes 
on UV Ceti stars. They conclude that less than 10 % of the ^Li observed in the interstellar 
medium are caused by stellar flares, de la Reza et ah (1981) determined the Li abundance 
in flare stars. They did not hud a general relation between Li abundance and chromospheric 
activity. On the other hand, Tatischeff & Thibaud (2007) show that the ®Li could be produced 
by the reaction '^He(^He, p)®Li in situ by solar-like flares. Montes & Ramsey (1998) reported 
the Li line enhancement during a long-duration flare in active binary. They suggest this Li 
enhancement is caused by spallation reactions during the flare. However, there is no clear 
evidence of productions of Li by stellar flares. In addition, Giampapa (1984) demonstrated 
that the spots and plages could alter the observed Li line on the basis of solar observations. 
Although many attempts have been made to study the possibilities of Li production in stellar 
flares, it is still controversial. 

Here we show the lithium abundance in superflare stars and discuss it’s relation with 
stellar activity, rotation velocity, and age of stars. In paper I and H, we determined the 
atmospheric parameters of superflare stars with high dispersion spectroscopy and discuss the 
spots and activity for understanding the nature of superflare stars. This is the third paper of 
our spectroscopic studies of superflare stars. In Sect. 2 and 3, we describe the observations 
and analysis of Li abundance in superflare stars, and in Sect. 4 we show the behaviors of Li 
abundance as a function of stellar parameters, discuss the reasons of Li distribution and the 
possibility of Li production by superflares. In Sect. 5 we summarize the paper. 

2. Observations 

High dispersion spectroscopy was carried out for 50 solar-type superflare stars (Notsu 
et ah 2013a, 2015a), which were found by Maehara et ah (2012) and Shibayama et ah (2013), 
on 2011 Aug. 3 (S11B-137S), 2012 Aug. 6-8, Sep. 22-25 (S12B-111N) and 2013 June 23-24 
(S13A-045N) using Subaru High Dispersion Spectrograph (HDS ; Noguchi et ah 2002). We 
applied the image sheer (jj2) in the observations of S13A (Tajitsu et ah 2012). In addition, we 
have observed 10 bright solar-type stars and the Moon for comparison. Among them, 8 stars 
are reported as “solar-twin” stars by the previous studies (King et al. 2005, Takeda &: Tajitsu 
2009, Datson et ah 2012). The details of observations and target stars are shown in Paper 1. 

The spectrum covers 6100-8820 A with a resolving power {R = A/AA) of 97,000 for 
SllB, 51,000 for S12B and 80,000 for S13A by 2 x 2 on-chip binning. This spectral range 
includes lines of Li I (6708 A), Ha (6563 A), and Ca H triplet (8498, 8542, 8662 A). The 
reduction was carried out in a standard manner using the IRAF echelle package^. We have 
excluded the binary stars in the following analysis. The 16 of 50 targets were recognized to be 

^ IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science 
Foundation. 
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binary systems which show double (or multiple) lines, variations of radial velocity or having 
visual companion stars (cf., Paper I). 

3. Analysis 

We derived the Li abundances and nsini (projected rotational velocity) from those 
spectra using the atmospheric parameters determined by Paper 1. The method of estimations 
of vsini is described in Notsu et ah (2013a) and Paper I. The objects are distributed in the 
range of Teg from 5,000K to 6,300K and logg from 3.5 to 4.9. We note that there are large 
differences between atmospheric parameters taken from Kepler Input Catalog (KIC : Brown et 
al. 2011) and our obtained values (Paper I). The estimated values of atmospheric parameters 
by spectroscopy are more reliable than that from photometric data such as KIC values (cf., 
Paper I). The adopted values in this analysis are shown in Table 1. 

For the abundance analysis, we used the analysis program SPTOOL^ which was devel¬ 
oped by Y. Takeda (private communication), based on Kurucz’s ATLAS9/WIDTH9 (Kurucz 
1993). We assumed local thermodynamic equilibrium (LTE) and derived abundances using 
the synthesis spectrum with interpolated model atmospheres taken from Kurucz (1993). The 
line data around Li I 6708 A region are taken from the list of Takeda & Kawanomoto (2005), 
which are based on the list of Smith et al. (1998) and Lambert et al. (1993). We assumed the 
contributions of ®Li are negligible. Figure 1 shows the spectra of all sample stars in the portion 
of the Li I 6708 line region. We used the value of solar abundances obtained by Asplund et al. 
(2009). The upper limit of Li abundance is estimated by the method of Takeda & Kawanomoto 
(2005). They used the averaged FWHM (full width at half maximum) of weak Fe lines and 
signal to noise ratios (S/N) for deriving the upper limit value. Adopted atmospheric parameters 
and derived Li abundances are shown in Table 1. 

We assume that there are two main sources of error. One is a systematic error by 
atmospheric parameters, and another is a random error depending on S/N. Systematic errors 
are caused by the choice of the model atmosphere. The typical estimated errors for Teg, log^f, 
Uturb, and [Fe/H] are 50 K, 0.2 dex, 0.2 km s“^, and 0.1 dex on the basis of Paper 1. Changes 
in the abundance of Li caused by errors of the above parameters are 0.05, 0.01, 0.06, and 0.00 
dex for the typical atmospheric parameter. The total systematic error is 0.08 dex, which is 
estimated from the root sum square of the above 4 systematic error values. 

The uncertainties arising from prohle htting error depend on the S/N and equivalent 
width of Li I line. We estimate the errors of measurement of equivalent width by using the 
formula derived by Cayrel (1988), which is taking in the S/N, resolution, and FWHM of line. 
Then, we estimate the severity of the value of A (Li) from that estimation. In most cases, the 
changes of A(Li) are smaller than 0.05 dex. However, a few objects show larger than 0.1 dex. 

^ http://optik2.mtk.nao.ac.jp/~takeda/sptool/ 
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The largest value is 0.3 dex for KIC 7354508. It is necessary to keep in mind that the low A(Li) 
stars have larger error than that of high ^(Li) stars. 

In the case of typical atmospheric parameter and Li abundance (Tes = 5,500K, log^f = 
4.4, ^(Li) = 2.1) in our sample, the effect of non-LTE is about 0.06 dex by using the grid of 
corrections (Carlsson et ah 1994). However, it must be noted that the correction of-0.3 dex is 
the most effective case for KIC 11610797, which has the highest and Li abundance in our 
sample. 

We then assume the typical error of Li abundance is about 0.15 dex in those objects 
from what has been discussed above. 

4. Results and Discussion 

4 ..I. Behaviors of Li abundance with effective temperature (Tes) 

Figure 2 shows the behavior of the Li abundances as a function of T^g of our target stars 
with F,G,K type stars (Takeda & Kawanomoto 2005). The Li depletion is remarkably seen 
in the stars whose temperature is lower than the Sun (Tgff ^ 5,500K). On the other hand, the 
main sequence stars having a temperature higher than 6,000K shows no Li depletion. This is 
because the evolution of the convection zone in the stellar atmosphere make transport the Li 
to deeper zone, and destruction of Li in the lower T^g stars. The depletion of Li in the stellar 
surface by convective mixing increases with a lapse of time. Hence, we can consider that high 
Li stars are young stars. 

However, a large diversity (by more than 2 dex) of the Li abundance is seen in stars with 
the solar temperature, despite the similarity of stellar parameters (cL, Takeda & Kawanomoto 
2005). Especially, the solar Li abundance is quite low. The standard stellar model cannot 
explain this diversity. The behavior of the Li abundance in solar-analog stars is still unclear. 

The Li abundances of superflare stars do not show clear relations with T^g. The temper¬ 
ature distribution of our target stars is slightly lower than that of the Sun. We have selected the 
target of superflare stars having the parameters similar to the Sun based on the temperature of 
KIC. The estimated temperature by KIC is systematically about 200K lower than other studies 
(e.g., Pinsonneault et al. 2012, Paper I). 

We can hnd the trend of low Li for rotation period (P) > 10 days and high Li for P < 
10 days. This trend is the same for the superflare stars having the parameters similar to the 
Sun (Figure 2b). The stellar rotation also might affect the Li depletion in superflare stars. 

We can consider that three stars (KIC 11610797, KIC 9652680, and KIC 8429280) which 
show an especially high value of Li are very young stars. It is also supported that these stars 
are young, because they show large projected rotational velocity (usini) and small rotation 
period (P). In addition to these three stars, about half of the objects show high Li compared 
with the stars in the Hyades cluster (Figure 3). The estimated age of Hyades cluster is 6.25 x 
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10® yr (e.g., Perryman et al. 1998), which is one of the old open clusters but younger than the 
solar age. It is acceptable to suppose that the young stars show high activities and superflares. 

However, more than ten stars with superflares do not show high values of Li (including 
upper limits) compared with the Hyades cluster. In particular, some objects show quite low 
values of Li like the Sun. Among them, stellar parameters of KIC 9766237 and KIC 9944137 
are also very similar to the Sun (Nogami et al. 2014). 

4-2. Rotational velocity and Li abundance 

Maehara et al. (2012) found the frequency of superflares increases with increasing the 
stellar rotation velocity. In general, the stellar rotation has a correlation with the age and activ¬ 
ity of the star. The young stars have rapid rotation, high activity and high Li abundance, but 
they have the decrease in the rotation, the activity and Li abundance with the age (Skumanich 
1972). 

Takeda et al. (2010) proposed that the stellar rotation (usini) may be the most important 
parameter in determining the surface Li content in the solar-analog stars. Figure 4 shows 
the correlation of Li abundance with usini of superflare stars and ordinary solar-analog stars 
(Takeda et al. 2010). The usini of solar-analog stars distributed mainly between from 2 to 5 
km s“^, and which show a clear correlation with Li abundance. 

Superflare stars show a wider range of usini than the solar-analog stars. We must note 
that this is not the typical distribution of superflare stars due to the selection bias (Paper I). 
However, we can hnd the trend of low Li and large usini in superflare stars compared with 
solar-analog stars. Some of them have a lower temperature (Tefj < 5,500K) and gravity (log^f 
< 4) than those of the Sun. Those stars tend to show lower Li abundances since they have the 
convection layer deeper than the typical solar-analog stars (5,500 K < Teg < 6,000K, log^f > 
4). It should also be added that those stars should be rapid rotation stars (usini > 5 km/s and 
P < 10 days). 

Exception is KIC 11764567, which shows large usini and low Li, but it has Sun-like 
atmospheric parameters. The rotational velocity of this star estimated from P is inconsistent 
with usini (cL, Paper H). This star might belong to a binary system. In the case of a binary 
system, observed superflares do not always reflect the primary solar-type star’s phenomenon. 

The Li abundances of KIC 4831454 and KIC 7354508 are remarkably deviated from the 
solar-analog stars. KIC 4831454 has a small inclination angle (Paper H), which should be rapid 
rotation by the estimation from P. On the other hand, KIC 7354508 shows lower Li abundance 
than the solar one in spite of slightly larger vsini than the solar one. We also found that the 
star shows a large value of radial velocity ( -112.8 km s“^ ; Paper I ). It is quite likely that this 
star is also belongs to a binary system. 

We can hnd that there are some superhare stars with low Li abundance (A(Li) < 1.5) 
and small usini ( ^ 2 km s“^) like the Sun. Among them, KIC 1197517, KIC 8359398, 
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KIC 11303472, and KIC 6504503 are not Sun-like stars. KIC 1197517, KIC 8359398, and KIC 
11303472 have low temperature (Teg < 5,300K) and KIC 6504503 has low gravity (log^f = 
3.6). However, there are some Sun-like superflare stars which show slow-rotation and low Li 
abundance like the Sun. 

In Paper I, we tried to estimate the age of superflare stars from isochrone. It is difficult 
to distinguish whether slightly low logg stars are post-MS or pre-MS only from the isochrone. 
Li abundance gives important clues about the age of the star. The Li abundances in many 
superflare stars have higher values than that of the Sun. Especially, those stars with a very high 
Li abundance are suggested to have a rapid rotation velocity from large nsini and brightness 
variations having a short period. The stellar rotation is also an index of the age. It seems 
reasonable to assume that the age of those stars are young. However, some superflare stars 
show small nsini (slow rotation) and low Li abundance. This result may indicate that those 
stars are not young. Superflare stars are not necessarily young on the basis of our spectroscopic 
observations. 

4-3. Ca II (ro(8542)) indexes and Li abundances. 

The ro(8542) index is residual core flux normalized by the continuum level at the line 
core of the Ca H (8542 A). The ro(8542) index reflects the activity of stellar chromosphere (e.g., 
Linsky et ah 1979) and well correlates with the intensity of the stellar mean magnetic field (cL, 
Paper H). Takeda et ah (2010) shows the positive correlation among Li abundances, ro(8542) 
and nsini for solar-analog stars. They conclude the depletion of surface Li in solar-type stars 
operates more efficiently as stellar rotation decelerates. 

In Paper H, the measured ro(8542) index shows that the superflare stars have higher 
chromospheric activity compared with the Sun. In addition to this, we found the correlation 
between the amplitude of the brightness variation and the ro(8542). This correlation can be 
explained by the difference of the spot size of superflare stars. 

Figure 5 shows the Li abundances as a function of ro(8542) index for solar-analog stars 
(Takeda et ah 2010). We could not find the positive correlation between Li abundance and 
ro(8542) of superflare stars compared with solar-analog stars. Duncan (1981) investigated the 
Li abundances and core flux of Ca H K line in solar-type stars. They found there are a number 
of high Li abundances and very little chromospheric emission flux, and the converse is rare. 
However, there are some low Li and high ro(8542) stars, which are not rare in superflare stars. 

These results indicate that some superflare stars have large active region but low Li 
abundance. Our results are not consistent with the Skumanich’s law (Skumanich 1972). The 
Skumanich’s law cannot apply to these superflare stars. We cloud not determine whether those 
stars are young or old from Ca H flux. 

The large spots are constituted to make enhancement of Li abundance than the normal 
region. Pallavicini et al. (1987) suggested that the strong Li line could be due to large cool 


spots in RS CVn stars. Superflare stars should have large spots on the stellar surface estimated 
from the amplitude of brightness variations. However, we can not hnd the strong Li line in 
such superflare stars with large spots. That effect may not be seen on about region of 10% 
(e.g., Notsu et ah 2013b) in the stellar surface. 

%% The possibility of Nucleosynthesis of Li in a Superflare. 

Tatischeff & Thibaud (2007) showed the possibilities of Li (®Li) production by stellar 
flares. However, there is no conclusive evidence that Li is produced in situ by stellar flares. 
Superflare stars are good objects to investigate whether flares can make Li or not. If Li produc¬ 
tion occurs during superflares, we could see an increase in the Li abundance of the superflare 
stars with frequent large flares. Contrary to such hypothesis, the results indicated that the Li 
abundance decreases as the number of superflares increases. Figure 6 shows the Li abundance 
as a function of the number of superflares and frequency of superflares in our target stars. We 
must note that the numbers of superflares should be a lower limit. It is because the detection 
of relatively small superflares is somewhat affected by stellar brightness variations (Shibayama 
et al. 2013). The frequency of superflares is obtained by dividing the number of superflares by 
the total observation time. KIC 8547383, KIC 7420545, KIC 6934317, and KIC 11764567 had 
superflares many times (> 30) in spite of their not high value of Li abundances (%(Li) < 2). 

The maximum energy of the superflare also does not show any correlation with the Li 
abundance. For example, KIC 10471412 shows the largest maximum energy of the superflare 
(5.2 X 10^^ erg : Paper I), but Li abundance is not high (%(Li) < 1.5). 

In our observations, we could not hnd any evidence of Li production by superhares. 
However, our observations could not show the small contributions of ®Li. Further research on 
Li isotope abundances of superhare stars would clarify the Li production by stellar hares. In 
order to know the reason for the high Li abundance of young stars, it is important to investigate 
the ratio of ®Li to ^Li in the solar-analog superhare stars. 

5. Summary 

We have estimated the Li abundance of superhare stars and investigate the correlations of 
Li abundance with stellar parameters, and the possibility of Li productions in stellar hares. Our 
spectroscopic observations show the slightly young solar-type stars tend to produce superhares, 
but old superhare stars exist. There is a possibility that superhares would be generated on 
our Sun. We could not hnd any evidence of nucleosynthesis of Li in stellar hares from our 
observations. 
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Fig. 1. Spectra around Li I 6708 line of the 34 superflare stars that show no evidence of binarity, 10 
comparison stars, and the Moon. The wavelength scale is adjusted to the laboratory frame. Co-added 
spectra are used here in case that the star was observed multiple times. 
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Fig. 2. Lithium abundances (A(Li)) vs. effective temperatures (Teff). a : Filled circles and squares 
indicate the Li abundance of superflare stars with P < 10 and P > 10 days respectively, and triangles also 
indicate the upper limit of Li abundance of superflare stars. Open small circles indicate the ordinary F, 
G, K-type main sequence stars (Takeda & Kawanomoto 2005), and open large circles indicate comparison 
stars in this observation, b ;A11 symbols are the same as a, but the data points of superflare stars are 
limited to solar-analog (5,500K < Teff < 6,000K, \ogg > 4) superflare stars. 
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Fig. 3. A(Li) vs. Teff with stars in the Hyades cluster. Upper panel : The symbols of filled circles and 
squares are the same as in Figure 1, but crosses indicate the stars of the Hyades cluster (Takeda et al. 
2013). The age of Hyades is 6.25 x 10® yr (e.g., Perryman et al. 1998). 
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Fig. 4. A(Li) vs. 'usinz. a : The symbols of filled circles, triangles, squares, and diamond indicate the 
measurements of this work. The plot of small open circles are taken from Takeda et al. (2010), which 
are the solar-analog stars. The small inclination angle stars (z < 20 ; Paper II) are shown as squares 
and diamond, b : All symbols are the same as a, but the data points of superfiare stars are limited to 
solar-analog (5,500K < T^s < 6,000K, logg > 4) superflare stars. 
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Fig. 6. a ; vl(Li) vs. the number of superflares (Shibayama et al. 2013). The symbols of filled circles and 
triangles indicate the measurements of this work. The triangles indicate the upper limit of Li abundance 
of superflare stars, b : A(Li) vs. the frequency of superflares (the number of superflares divided by total 
observation time). 
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Table 1. Stellar parameters, estimated rotation period, r^sini, Li abundance, ro index, number of superfiares, and total obser¬ 
vation time. 


KIC ID 

Teffl 

[K] 

log 5^ 

[cm s“^] 

Vt^ 

[km s“^] 

[Fe/H]i 

7l(Li) 

[km s“^] 

pi 

[days] 

ro (8542)2 

N_f® 

T.obs^ 

[days] 

Remarks 

3626094 

6026 

4.15 

1.17 

-0.03 

2.56 

2.9 

0.7 

0.24 

6 

483 

i < 20® 

4742436 

6033 

4.20 

1.09 

-0.15 

2.56 

3.6 

2.3 

0.30 

14 

483 

i < 20® 

4831454 

5672 

4.59 

1.13 

0.04 

2.83 

2.1 

5.2 

0.50 

4 

483 

i<20 

6503434 

5902 

3.63 

0.96 

-0.20 

2.37 

5.5 

3.9 

0.27 

13 

389 


6504503 

5484 

4.56 

0.77 

0.23 

<1.2 

1.6 

31.8 

0.25 

1 

389 

i < 20® 

6865484 

5800 

4.52 

1.11 

-0.12 

2.33 

2.7 

10.3 

0.48 

12 

473 


6934317 

5694® 

4.42® 

0.87® 

-0.03® 

1.25® 

1.9® 

2.5® 

0.33® 

48 

483 

i < 20® 

7093547 

5497 

4.62 

0.18 

0.26 

1.04 

2.1 

14.2 

- 

1 

473 


7354508 

5620 

4.09 

0.92 

-0.10 

0.81 

3.2 

16.8 

0.40 

1 

473 


7420545 

5355 

3.56 

1.30 

-0.03 

1.57 

6.3 

36.2 

0.40 

33 

483 


8359398 

5214 

4.85 

0.75 

-0.13 

1.11 

1.8 

12.7 

- 

3 

473 


8429280 

5055^ 

4.41^ 

1.00^ 

-0.02^ 

2.9^ 

37.1^ 

1.2 

0.68 

6 

483 


8547383 

5606 

4.38 

0.83 

-0.02 

<2.0 

4.4 

14.8 

- 

40 

473 


8802340 

5529 

4.43 

1.25 

-0.13 

1.92 

4.3 

10.3 

0.49 

2 

483 


9412514 

6075 

3.72 

1.55 

0.10 

1.58 

8.1 

3.7 

0.22 

1 

483 

i < 20® 

9459362 

5458 

3.82 

1.39 

-0.78 

1.63 

7.0 

12.6 

- 

19 

473 


9583493 

5605 

4.34 

1.28 

-0.12 

2.61 

7.2 

5.5 

0.51 

1 

483 


9652680 

5618 

4.80 

1.00 

-0.30 

3.39 

34.2 

1.5 

0.68 

2 

483 


9766237 

5606® 

4.25® 

0.88® 

-0.16® 

<1.4® 

2.1® 

14.2 

0.26® 

1 

473 


9944137 

5666® 

4.46® 

0.93® 

-0.10® 

<1.4® 

1.9® 

12.6 

0.27® 

1 

473 


10252382 

5334 

4.01 

1.40 

-0.61 

2.04 

5.5 

16.8 

- 

3 

473 


10387363 

5381 

4.75 

0.65 

0.07 

<1.5 

1.5 

29.9 

- 

1 

473 


10471412 

5776 

4.53 

0.91 

0.15 

<1.4 

2.7 

15.2 

0.35 

4 

483 


10528093 

5180 

3.96 

0.71 

-0.14 

1.84 

11.4 

12.2 

0.67 

18 

483 


11140181 

5552 

4.60 

1.25 

-0.09 

1.49 

3.6 

11.5 

0.56 

1 

483 


11197517 

5284 

4.59 

0.62 

0.42 

<1.0 

0.9 

14.3 

0.19 

1 

483 


11303472 

5193 

4.57 

0.93 

-0.14 

<1.0 

2.1 

13.5 

0.49 

4 

473 


11390058 

5921 

4.43 

1.11 

-0.15 

2.50 

3.5 

12.0 

0.32 

3 

418 


11455711 

5631 

3.85 

1.16 

-0.23 

2.15 

7.9 

13.9 

0.40 

3 

473 


11494048 

5907 

4.49 

1.03 

0.01 

2.46 

3.4 

14.8 

0.32 

1 

408 


11610797 

6209 

4.41 

1.70 

0.26 

3.61 

22.5 

1.6 

0.63 

4 

483 


11764567 

5736 

4.02 

1.71 

-0.08 

1.69 

14.7 

22.4 

0.49 

34 

483 


11818740 

5500 

4.84 

0.80 

0.00 

<1.9 

3.0 

15.4 

- 

1 

473 


12266582 

1 rpi_X-T 

5576 

4.61 

1.20 

-0.06 

2.34 

4.8 

6.9 

0.52 

1 

483 



^ The atmospheric parameters and brightness variation periods (P) which are estimated in Notsu et al. (2015a : Paper I). 
^ Normalized intensity of the line center of Ca II 8542 line takL§ from Notsu et al. (2015b : Paper II). 

^ N_f is the number of superflares (Shibayama et al. 2013). 


T_obs is the observation time of each superflare stars by Kepler (QO - Q6). 
^ Notsu et al. (2015b :Paper II). 


Notsu et al. (2013b). 





Table 2. Stellar parameters, estimated rotation period, r^sinz, and Li abundance. 


ID 

[K] 

loggi 
[cm s“^] 

[km s“^] 

[Fe/H]i 

[km s“^] 

A(Li) 

59 Vir 

6127 

4.32 

1.38 

0.19 

6.2 

2.92 

61 Vir 

5581 

4.52 

0.89 

0.01 

1.2 

<3.0 

18 Sco 

5760 

4.37 

1.02 

0.03 

2.3 

1.50 

HD163441 

5768 

4.45 

0.87 

0.05 

2.5 

1.55 

HD173071 

5982 

4.41 

0.98 

0.18 

2.7 

2.48 

HIP71813 

5786 

4.30 

1.07 

0.01 

2.3 

0.49 

HIP76114 

5746 

4.50 

0.98 

0.02 

1.8 

0.76 

HIP77718 

5830 

4.37 

0.97 

-0.11 

2.5 

2.28 

HIP78399 

5830 

4.44 

1.04 

0.07 

2.5 

1.83 

HIP100963 

5834 

4.56 

1.07 

0.01 

2.4 

1.70 

Sun(Moon) 

5783 

4.44 

0.85 

0.03 

2.4 

< 1.0, 0.92^ 


j. iifcj SLfcJiicii pciiciiiitJiitJia ciit; Lajs.fc;ii iiuiii inulhu et a,i. i^^ux<jci.jra,pei l ). 

^ The Li abundance of solar value are taken from Asplund et al. (2009). 
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